Sequential Landau-Zener Transitions in A Spin-Phonon Coupled System 
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We investigate the Landau-Zener process of a harmonically trapped spin-i particle with spin- 
phonon coupling. Because of the spin-phonon coupling, the energy level structures are modified by 
the Franck-Condon effects, in which some avoided energy-level-crossings are almost closed and some 
energy- level-crossings are opened. The close of avoided energy-level-crossings and the open of energy- 
level-crossings result in the Franck-Condon blockade and the vibrational transitions, respectively. If 
the spin-phonon coupling is sufficiently strong, the sequential Landau-Zener transitions of ladder- 
like population transition appear in the Landau-Zener process, in which an initially ground-state 
system is driven from one far-off-resonance limit to the other far-off-resonance limit. We derive a 
generalized Landau-Zener formula which is well consistent with our numerical results. 
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Landau-Zener (LZ) problem [l], 0], a well-known fun- 
damental problem in time-dependent quantum mechan- 
ics, concerns how non-adiabatic transition appears in 
a two-level system driven through an avoided energy- 
level-crossing. According to the quantum adiabatic the- 
orem @, 3, if the system varies infinitely slow, non- 
adiabatic transition will not take place and the system 
will always be in an eigenstate of its instantaneous Hamil- 
tonian. The studies of LZ transition are not only of great 
fundamental interests but a l so °f extensive appli- 

cations I3"17| m quantum state engineering, quantum 
interferometry and quantum computation etc. 

Up to now, there are lots of theoretical and experimen- 
tal studies of LZ transitions in systems of decoupled in- 
ternal spin states and external motional states. However, 
to the best of our knowledge, there are almost no stud- 
ies on LZ transitions in closed quantum systems whose 
internal spin states and external motional states are cou- 
pled. Spin-phonon coupling (SPC) is a typical coupling 
between internal spin states and external motional states. 
How SPC affects a LZ process ? What does SPC bring to 
a conventional LZ transition? 

It has been demonstrated that strong SPC may in- 
duce Franck-Condon (FC) blockade and vibrational side- 
bands. The FC principle [H HI] tells us that the FC 
blockade takes place if the FC factor, which is defined 
as the square of the overlap integral between the vi- 
brational wave-functions of the two involved states, is 
sufficiently small to be ignorable. The FC blockade in- 
duced by strong SPC has been found in several systems, 
such as molecular junctions [20l-l22l| . nano-tube quantum 
dots neutral atoms [23], and a single trapped 



ion [28(. On the other hand, nonzero FC factors may 
cause vibrational sidebands 0, El H3, [H-Hl] , that 
is, the population transfer or the electronic tunnelling 
has been shown to excite vibrational modes. Is there any 



signature of the FC effects in the LZ process of a spin- 
phonon coupled system? 

In this Letter, we investigate the LZ process of a spin- 
phonon coupled system of a single spin-^ particle in 
a spin-dependent harmonic potential. We explore how 
SPC affects energy-level structures and the LZ process. 
The energy gaps of avoided energy- level-crossings become 
narrow when the SPC becomes strong. At the same 
time, some energy-level-crossings are gradually opened 
and then closed. The appearance of FC blockade and 
vibrational sidebands are direct results of the close of 
avoided energy-level-crossings and the open of energy- 
level-crossings, respectively. Under sufficiently strong 
SPCs, in contrast to the LZ transition in a system with- 
out SPC, the sequential LZ transitions of ladder-like pop- 
ulation transition appear. We find that the FC blockade 
corresponds to the absence of some specific population 
steps. By treating the sequential LZ transitions as a se- 
quence of conventional two-level LZ transitions and ap- 
plying the conventional two-level LZ formula again and 
again, we obtain a generalized LZ formula for the final 
populations. Based upon the experimental techniques 
for exploring the FC blockade, it is possible to test our 
prediction of the sequential LZ transitions in future ex- 
periments. 

We consider a harmonically trapped particle, which 
may be an electron, an atom, a molecular or a quantum 
dot and only two of its internal spin states of different 
magnetic dipole moments are involved. Below, we call 
this particle as a spin-i particle and denote two involved 
spin states as spin-up state |f) and spin-down state 
In the LZ process, the two spin states are coupled by 
lasers with a linearly sweeping detuning. For simplicity 
and without loss of generality, we concentrate our studies 
on one-dimension systems with SPCs induced by gradi- 
ent Zeeman fields. The Hamiltonian for our spin-phonon 
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coupled system reads as, 

H = H LZ + H ho + H s 



fective Hamiltonian, 



where H LZ = -^-a z - ^fa x with 5{t) 
conventional LZ Hamiltonian, H^o - 



5q + at is the 



2m oz 2 2 £ 

describes the external motion and H spc — mui 2 Zoza z 
characterizes the SPC. Here, a XiZ are Pauli matrices, 
h is the Planck constant, m is the particle mass, Q is 
the Rabi frequency, S(t) is the detuning, and zq is the 
SPC strength. Taking into account the SPC, the particle 
feels a spin-dependent potential V<j{z,t) = -~^moj 2 z 2 

with (a z ) 



l_,2^ (CTz) , o) 2_M£) ((T 
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and (<j z ) = — 1 for a =].. In Fig. 1, we show the schematic 
diagram for the LZ process of Hamiltonian (1). 







FIG. 1: Schematic diagram for the Landau-Zener process of 
a spin-i particle in a harmonic trap. The red dashed curves 
and the blue solid curves denote the spin-dependent harmonic 
trap Vf(z) and Vj_(z), respectively. The detuning is linearly 
swept from one far-off-resonance limit (left column) to the 
other far-off-resonance limit (right column) through the zero 
point (middle column), (a) Zo = for the case of no spin- 
phonon coupling and (b) zq 7^ for the case of spin-phonon 
coupling. 



Denoting the quantum state of Hamiltonian (1) as 
\^>) = ( ^(^'t) J ' ^ S time-evolution obeys the coupled 
Schrodinger equations, 





' h 2 d 2 




2m dz 2 




' h 2 d 2 


dt 1 


2m dz 2 



+ V t (z,t) 
+ V i (z,t) 



-0 t - 



Ml 
~2~ 

m 



it, 



(2) 
(3) 



with J dz (\ip-f(z, t)\ 2 + \4>\.{z, t)\ 2 ^j = 1 imposed by the 

normalization condition. Given the eigenstatcs \un) 

for V a (z,t), the wave- functions can be expressed as 

+00 

^a{z,t) = J2 C an {t)\an) with C an {t) = {na\ijj a (z,t)) . 
In the basis composed of \cm), the system obeys an ef- 



(!) H = J^Ean Wn) (na\ - ^E^ t [|| n') (n f| + h.c] 



cr,n n\n 

(4) 

Here, the FC factors F n , n = \(n' || f n)\ 2 = \(ri t| I n)\ 2 
and E an = nhu z — (a z )hS/2 with (a z ) = +1 for a =f 
and (a z ) = — 1 for a =\.. Comparing with Hamiltonian 
(1), the zero-energy point of Hamiltonian (4) is shifted to 
(ftw z — mui 2 z 2 ^ /2. From the effective Hamiltonian (4), 
the complex amplitudes C an (t) obey 

ih^- t C tn (t) = E tn C tn (t) - ^^V^CWi), (5) 

n' 

rt§;C in (t) = E in C in (t) - ^^/F^C t „(t). (6) 

n' 

Obviously, if F n > n — 0, there will be no population trans- 
fer between \an) and \<r'n') with a ^ a'. 

Below, we analyze the time evolutions and the pop- 
ulation transitions in the LZ process. By numerically 
integrating the coupled Schrodinger equations (2) and 
(3), we obtain ip a {z,t) and then calculate both the spin 
populations P a {t) = J dz \ipa-{z, t)\ and the phonon pop- 
ulations P an (t) — \C an (t)\ 2 . In our numerical simulation, 
we have chosen the natural units of m = 1, h = 1 and 
u> z = 1- The initial state is chosen as the ground state in 
the negative far-off-resonance limit (S < and \S\ ^> 17), 
that is, ipi(z,t = 0) = || 0) and ip^(z,t = 0) = 0. In 
the LZ processes, the Rabi frequency SI is fixed as 0.2, 
the detuning is linearly swept according to S = Sq + at 
with the initial detuning <5o = — 1 and the sweeping rate 
a = 2.5 x 1(T 4 . 
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FIG. 2: Time evolution of probability densities in the Landau- 



Zener process. The 
column is |^i-|-(2,f)| 2 . 



left column is |^(z,t)| and the right 
The top, middle and bottom rows cor- 
respond to the SPC strength zo = 0, 1.5 and 2.5, respec- 
tively. The other physical parameters are chosen as So = —10, 
a = 2.5 x 1(T 4 and O = 0.2. 



In Fig. 2, for different values of the SPC strength zq, 
we show the time evolution of the probability densities 
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|V't,i( z ! ^)l in the LZ process. If there is no SPC, the 
system undergoes adiabatic evolution and the two prob- 
ability densities keep in Gaussian shapes. The unchanged 
shapes of \ip^ t ±(z, t)\ 2 for z Q = indicate the absence of 
phonon excitations. When the SPC strength zo increases, 
although \ipi(z,t)\ 2 keeps in a Gaussian shape, multi- 
hump structures gradually appear in t)\ . The ap- 
pearance of multi-hump structures in |i/^(z,i)| 2 is a sig- 
nature of phonon excitations. In particular, the signifi- 
cant changes of the two probability densities sequentially 
take place in the vicinity of S = nftiu (where n are non- 
negative integers). 
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FIG. 3: Sequential Landau-Zener transitions induced by spin- 
phonon coupling, (a) The populations Pj, (red curves) and 
Pf (blue curves) versus the detuning S for zo = and (b) 
the ones for zq = 2.5. (c) and (d) are the two energy spectra 
corresponding to (a) and (b), respectively. The sequential 
Landau-Zener transitions in (b) have ladder-like population 
transition with steps at S = nhw z . The numbers 1 to 7 in (d) 
label the avoided energy-level-crossings, which correspond to 
the population steps in (b). 

In Fig. 3 (a) , we show the population transitions corre- 
sponding to the time evolutions in the top row of Fig. 2. 
Due to the absence of SPC, the phonon states are spin- 
independent and the FC factors F n ' >n are non-zero if and 
only if n' = n. Therefore, the time evolution of spin 
states and phonon states are decoupled and phonon ex- 
citations will not take place in the LZ process. In the 
corresponding energy spectrum, the avoided energy-level- 
crossings only appear around (5 = and these avoided 
energy-level-crossings dominate the population transfer 
in the LZ process, see Fig. 3 (c). In the LZ process, 
as labeled by the arrows, the system evolves along its 
instantaneous ground state due to the sweeping rate 
a = 2.5 x 1CP 4 is sufficiently small. 

In Fig. 3 (b), we show the population transitions cor- 



responding to the time evolutions in the bottom row of 
Fig. 2. Due to the strong SPC, the FC blockade and 
the sequential LZ transitions appear in the LZ process. 
Corresponding to the significant density changes in Fig. 
2, a series of population steps appear at S = nhuj. This 
ladder-like population transition is a direct signature of 
the sequential LZ transitions. In addition, unlike the 
conventional LZ transition, there is no population step at 
5 = 0. The absence of the population step at S = is a re- 
sult of the FC blockade between the lowest-phonon states 
|4- 0) and |t 0). Based upon our numerical results for dif- 
ferent values of zq (not shown here), the FC blockade 
appear only when z Q is sufficiently large and more pop- 
ulation steps will disappear due to the FC blockades be- 
tween the lowest-phonon state |t 0) and the high-phonon 
states |t n) take place for larger z . In the energy spec- 
trum, because of the FC effects, some avoided energy- 
level-crossings are almost closed due to the correspond- 
ing FC factors are very small, and some energy-level- 
crossings opened due to the corresponding FC factors 
become non-zero, see Fig. 3 (d). Therefore, as labeled by 
the arrows, the system undergoes sequential LZ transi- 
tions in which sequential phonon excitations accompany 
the ladder-like spin population transition. 

To characterize the phonon excitations, one has to an- 
alyze the population dynamics of \an). Given ip a (z,t) 1 
it is easy to obtain the phonon populations P an {t) = 
\C an (t)\ 2 with C an (t) = {na\i]j a {z,t)). In the absence of 
SPC, the phonon populations P n (t) — P± n (t) + P-f n {t) 
keep unchanged, i.e. P n (t) = P„(0). In the presence 
of SPC, due to the occurrence of the phonon excitations 
from |t 0) to |t n) around 5 = nfkj z , step- like changes 
appear in the phonon populations. In other words, the 
population step at 5 — nhu z is caused by the sideband 
transition between |t 0) and |t n). 

The sequential LZ transitions can be treated as a 
sequence of conventional two-level LZ transitions. To 
calculate the transition probability of a two-level LZ 
transition in the sequence, in addition to the sweeping 
rate a, one has to know the energy gap for its avoided 
energy-level-crossing. The energy gap sensitively de- 
pends on the SPC strength z . By diagonalizing the ef- 
fective Hamiltonian (4), the energy gap for the avoided 
energy-level-crossing between |t n) and |t n') is given as 
A n / n = fl\/F n 'n- If there is no SPC, A n > n = Q6 n ' n with 
S n 'n denoting a Kronecker delta function. If the SPC is 
sufficiently weak, Aoo is still large enough and the LZ pro- 
cess is similar to the conventional LZ process of no SPC. 
The energy gaps A„„ become narrow when the SPC be- 
comes strong and at the same time the other energy gaps 
A„<„ of n 1 7^ n are gradually opened and then closed. If 
the SPC is sufficiently strong, the almost vanishing A„< n 
may induce the blockade of the vibrational transitions 
|t n) |t n') in the LZ process. 

Now, by using the conventional two-level LZ formula 
again and again, we derive the analytical formula for the 
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final populations. The sequential LZ transitions can be 
decomposed into a sequence of conventional two-level LZ 
transitions: || 0) ^> |f0), || 0) <S> |f 1), •■•,14-0) 

■••, see Fig. 3 (d). Therefore, by applying the 
conventional two-level LZ formula one by one, the final 
populations are given as 

P t o = l-P L z(A 00 ,a), 
P n = P L z(A 00 , a)(l - P LZ (A 01 ,a)), 
P t 2 = P L z ( A Q0 , a) P L z { A 01 , a) (I - P LZ (A Q2 , a)), 



P 



7i f —n— 1 



II p Lz{A w ,a) 



(l-P LZ (A 0n ,a)), (7) 



Pi 



In 



n ' —-\-oo 

I] P L z(A 0n ,,a), for n = 0, 

n'=0 

0, for n ^ 0, 



with Plz{Aq ti , a) = cxp ^— j given by the conven- 
tional two-level LZ formula. 



TABLE I: Comparison between numerical and analytical re- 
sults of final populations. The first column is the SPC 
strength zo, the second column is the FC factor, the third 
column is the final population obtained from numerical inte- 
gration, and the last column is the final population given by 
the generalized LZ formula (7). The rows of zo = and 2.5 
correspond to (a) and (b) of Fig. 3, respectively. 



SPC 


FC Factor 


Numerical 


Analytical 


zo = 


Foo = 1 


P t0 = 0.9998 


P t0 = 1.0000 



zo 



TO 



"to 



Foi = 4.7 x 10" 



Pfi = 0.0059 P n = 0.0058 



F 2 

F 03 

F 4 

Fos 
F oe 
For 



0.0003 
0.0012 
0.0038 
0.0095 
0.0197 
0.0353 



Pf3 



0.0367 
0.1383 



Pn 

Pf3 



0.0357 
0.1355 



P t4 = 0.3152 P t4 = 0.3117 
P t5 = 0.3535 P t s = 0.3556 



T-5 
Pf6 
Pf7 



0.1378 
0.0121 



t5 
Pf6 

Pn 



0.1423 
0.0128 



The final populations given by the generalized LZ for- 
mula (7) are well consistent with our numerical results. 
In Table I, for different values of the SPC strength Zq, we 
compare the final populations estimated by the general- 
ized LZ formula (7) with the corresponding ones obtained 
by numerical integration. For zq — 0, which corresponds 
to the numerical results shown in Fig. 3 (a), the phonon 
number n is exactly truncated at because there is no 
SPC. For zq = 2.5, which corresponds to the numerical 
results shown in Fig. 3 (b), the phonon number n can be 
approximately truncated at 7. In particular, the FC fac- 
tor Foo is in order of 10~ 6 and the final population P-fo is 
in order of 10~ 4 . As the initial state is |i 0), such a small 
Pto is a signature of the FC blockade between |i 0) and 
|t 0). For a stronger SPC, there may appear more FC 



blockades between |i 0) and |t n) with n up to a larger 
integer number. The numerical and analytical data of 
the final populations clearly show that the largest abso- 
lute difference is less than 0.005 and the largest relative 
difference is less than 5.5%. This means that the gener- 
alized LZ formula (7) is a very good estimation for the 
final populations. 

By using a single electron spin in suspended car- 
bon nanotube quantum dot |24j, a single ion in spin- 
dependent Paul trap [28j|, or ultracold atoms in spin- 
dependent optical lattices |27j |. it is possible to test our 
prediction of the sequential LZ transitions. Here, for two 
examples, we briefly discuss the experimental possibil- 
ity via a sin gle trapped ion [28| or noninteracting ultra- 
cold atoms [27j with tunable SPCs. In the case of a 
single trapped ion, the spin-dependent potential can be 
formed by imposing a gradient magnetic field on the Paul 
trap [33[ and the SPC strength is determined by the field 
gradient. In the case of noninteracting ultracold atoms, 
the spin-dependent optical lattices [3] can be created by 
two polarized counter-propagating laser beams and the 
SPC strength can be adjusted by tuning the polarization 
angle. In both two cases, the two internal spin states can 
be coupled by Raman laser beams and the detuning can 
be varied by modifying the laser frequency. 

In summary, we have studied the LZ problem of a har- 
monically trapped spin-^ particle with SPC. The intrin- 
sic mechanism of the sequential LZ transitions and the 
direct signatures of the FC effects in the LZ process have 
been explored. Due to the FC effects, the energy gap for 
the avoided energy-level-crossing between |i n) and |t n') 
is given as A n > n — Q\/F n ' n . The sequential LZ transi- 
tions, which have ladder-like population transition, take 
place when the SPC is sufficiently strong. The presence 
of population steps at S = nhuj z is a signature of the 
vibrational sidebands and the absence of some specific 
steps is a signature of the FC blockade. The sequential 
LZ transitions can be treated as a sequence of conven- 
tional two-level LZ transitions. We analytically derive a 
generalized LZ formula, which gives an excellent estima- 
tion for the final populations. We also briefly discuss the 
experimental possibility. 
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